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Abstract 

It is planned to use atomic processes to spread out most of the heating power over the first 
wall and side walls to reduce the heat loads on the plasma facing components in ITER to ~ 
50 MW. Calculations indicate that there will be 100 MW in bremstrahlung radiation from the 
plasma center, 50 MW of radiation from the plasma edge inside the separatrix and 100 MW 
of radiation from the scrape-off layer and divertor plasma, leaving 50 MW of power to be 
deposited on the divertor plates. The radiation losses are enhanced by the injection of 
impurities such as Neon or Argon at acceptably low levels (-0.1 % Argon, etc.) 

1 . Introduction 

Power exhaust will be a key issue for the next generation of fusion experiments such as 
ITErI . The alpha heating power in ITER will be about 300 MW, with possibly another 100 
MW of auxiliary power for driven operation. This gives a peak heat load of about 30 
MW/m 2 or larger. To maintain the surface temperature of the divertor plates at an acceptable 
level requires that the peak heat loads be reduced to ~ 5 MW/m 2 .which requires that the total 
heat loads on the divertor plates be 50 MW or less. This can be accomplished by using 
impurity radiation, enhanced by the injection of impurities such as Ne or Ar, to spread out 
the plasma energy on the first wall and divertor chamber walls. An analysis of the various 
channels for energy losses from the plasma for the plasma core, the plasma edge inside the 
separatrix (mantle), the scrape-off layer and divertor plasma, and the divertor plates 
indicates that it should be possible to radiate about 250 MW of power to the walls before the 
energy reaches the plate (Figure 1 and Table 1). 
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Table 1 . Power flow in ITER 



Bremsstrahlung 

Edge impurity radiation (inside separatrix) 
SOL/divertor impurity radiation (outside separatrix) + 
H and He radiation and charge exchange losses 
Power on target plates 


100 MW 
50 MW 
100 (up to 250) MW 

50 MW 


total heating power 


300 (up to 400) MW 



2. Detailed Power Losses in ITER 



Bremsstrahlung losses will be larger in ITER than in present experiments due to the higher 
electron density, temperature and energy confinement time. The ratio of the Bremsstrahlung 
emission P Brem to the alpha heating P a is independent of size and density (equation 1), where 

Cb is 2xl0 -40 Zeff MW/(m 6 eV - 5 ), n e and T e are the volume averaged density and 
temperature, fDT= n DT/ n e> E a = 3.5 MeV, and <ov>dt is the fusion reaction rate. For 
ITER, (<T>~10 keV, f a ~0.1, f Be ~0.01, f Ne ~0.003 and P B r em /Pa=0.33 (Fig.2). P B rem/Pa 
decreases with T e as Te" 1 - 5 , and increases with impurity level (eq.2) where fz=ni m purity/ n e- 
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Figure 1 . Power flow in ITER. 



Figure 2. Ratio of Bremsstrahlung radiation 
to alpha heating as a temperature for nne/n e 
= 0, 0.05, 0.1, 0.15 and 0.2. PBrem/Pheat is 
also plotted for P hea t = 20 MW for Dffl-D 
and 100 MW for ITER. 
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There will be impurity radiation from the main plasma edge (mantle) just inside the 
separatrix to the first wall due to line emission from intrinsic impurities and injected 
impurities such as Neon or Argon. This radiation can be calculated from the radial heat 
conduction equation and the radiation emission rates for candidate injected impurities 
(equation 3). If one assumes that the density and the electron transport coefficients are 
roughly constant, as is approximately true for H-mode plasmas ^, the radiation losses are 
proportional to an integral of the radiation emission rate^ '^(Figure 3). The effect of heat 
conduction on the radiation losses and the different scaling with density and impurity 
fraction is significant. With conduction, Qi ~ n Vf , which is much weaker than the 
volume loss rate ~ n f . 
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For the range of expected conditions for ITER for Zeff -1.6 with fHe = 0.1, fBe = 0.01, 



(3) 



8 x 10 19 m" 3 < n , < 10 20 m" 3 , 0.125 x 10 20 
edge 



nr 1 s" 1 < k < 2 x 10 20 nr 1 s" 1 , Neon, 



Argon and Krypton can provide radiation losses of about 50 MW from the plasma edge (Qj_ 
~ 0.05 MW/m 2 , Table 2). These analytic estimates are consistent with detailed calculations 
using the WHIST 1 1/2 tokamak transport code with a multi-species impurity transport 
package with the constraint that AP a / Pa ^ 5% andiE ~ 2 x xeITERP-89. 
Table 2 Radial Impurity Radiation losses for Be, C, Ne, Ar and Kr for ITER conditions 





Be 


c 


Ne 


Ar 


Kr 


Assumed fraction (Z e ff < 1.6) 


0.01 


0.01 


0.0033 


0.001 


.00025 


Q ± /(n e (K 1 f z )0-5)MWm* 


0.1 


0.4 


1.0 


3.0 


10 


Q ± (MW/m 2 ) for 5 x 10 iy < n e < 


0.002 to 


0.007 to 


0.01 to 


0.02 to 


0.03 to 


10 2 °and 0.125 <k< 2 


0.014 


0.06 


0.08 


0.13 


0.22 



*K ± in 10 20 nr 1 s" 1 , n e in 10 20 



nr 



A similar model can be used to estimate the radiation losses from the scrape-off layer and 
divertor plasma equation 4). The parallel heat conduction is classical and the pressure is 
constant along the field lines. The radiation efficiency can be computed as a function of the 
upstream temperature(Figure 5). The upstream temperature can be calculated by integrating 
the heat flux equation from the midplane to the radiation region. The impurity radiation 
losses can be computed for candidate ITER upstream densities (Table 3). For heating 
powers of 100 to 200 MW, Qn ~ 0.3 to 0.6 GW/m 2 for each divertor leg. Neon and higher 
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Z impurities can radiate 100 MW or more from the divertor and SOL plasma. These losses 
will be supplemented by H, He and Be radiation and charge exchange losses. These results 
are consistent with detailed 2 D B2 code divertor modelling calculations using a multi- 
species transport model for Neon or Argon impurities (Table 4). 

Power radiated from main plasma 
as a function of the peak Argon concentration 
WHIST (Mandrekas and Stacey) 





10 100 

T (eV) 

Figure 3. The normalized radiation losses 
y nJ o,y,=^2f o L z (T e )dT e vs.T e . 



0.06 0.08 0.1 0.12 0.14 

f z (%) 

Figure 4. Power radiated from the main 
plasma as a function of the peak Argon 

concentration from a Whist calculation. 
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Table 3 Divertor radiation efficiencies for DIII-D and ITER. 



Element 


Be C Ne Ar Kr 


f z (%) for Zeff ~ 1.6 with fRe = 0.1, fBe ~ 0.01 
QllITER/V(fz(%)/Z e ff) at 250 eV (GW/m 2 ) 
QllITER (coronal equihbrium) for n e ~ 10^0 m"3 


0.01 0.01 .0033 .001 .00025 
0.08 0.2 0.7 2 4 
0.065 0.16 0.33 0.52 0.52 



These radiation levels can be enhanced by charge exchange recombination and rapid 
recycling of the impurities^ the radiation losses calculated as before(Figure 6). To 
radiate the required 100 MW in the divertor plasma takes relatively modest levels of 
recycling and neutral densities (Table 5). 

The power balance in Table 1 thus appears to be achievable for operation with 300 
MW of ignited operation. Driven operation with 400 MW of total heating will require larger 
radiation levels from the divertor. 
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The conducted heat flux (recycling enhanced) vs T 
for n /n e = 10' 7 , 10 6 < n t T recyc < 10 14 s cm 3 




T (100 eV) 



Figure 5. Radiation integral for parallel 

. c , c . ■ , Figure 6. Recycling enhanced radiation 

transport tor a number of materials. b J b 



integral for parallel transport for carbon. 
Table 4 B2 calculations of the radiated power in the ITER divertor for two average Neon 

concentrations and upstream densities 



ne sep (10 19 m -3) 


fNe (%) 


Prad(MW) 


4 


0.5 


70 


4 


0.25 


40 


6 


0.5 


100 


6 


0.25 


60 



Table 5 Enhanced Be, C, Ne, and Ar divertor radiation efficiencies for DIII-D and ITER. 



Element 


Be C Ne Ar 


f z (%) for Zeff ~ 1.6 with fRe = 0.1, fee ~ 0.01 
QllITER/V(fz(%)/Z e ff) at 250 eV required to 
radiate 100 MW 

n Q /n e (required to radiate 100 MW) 
.3 

ne^recycle (s cm )(required to radiate 100 MW) 


1 1 0.33 0.1 
0.8 0.8 1.4 2.5 

6 x 10" 3 10" 3 3 x 10" 4 10" 3 
2xl0 10 8xl0 10 2x 10 11 10 10 
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